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Abstract 

In this paper, we present a collection of results on the observability of quantum mechan- 
ical systems, in the case the output is the result of a discrete nonselective measurement. 
By defining an effective observable, we extend previous results, on the Lie algebraic char- 
acterization of observable systems, to general measurements. Further results include 
the characterization of a 'best probe' (i.e. a minimally disturbing probe) in indirect 

| measurement and a study of the relation between disturbance and observability in this 

case. We also discuss how the observability properties of a quantum system relate to 

' the problem of state reconstruction. Extensions of the formalism to the case of selective 

£Nj , measurements are also given. 

jjy ', 1 Introduction 

The structural properties of controllability and observability have been studied in depth for 
. deterministic control systems (see e.g. [201) of the form 

> : 

x = f(t,x,u), (1) 

in ■ 

f^*) ■ with output 

as : 

; y = y{x). (2) 

In Q x is the state of the system varying on a given manifold M, u is the control, / 
a smooth vector field and y a smooth map M — > R which models how observations on the 
system depend on the state. For quantum systems, the study of controllability has received 
greater attention (see e.g. [2], |H], ^3], E|)- A study of the observability for quantum 
systems is complicated by the fact that, in general, the output has a probabilistic nature 
and the associated probability distribution depends on the current state. Moreover, different 
types of measurements can be considered according to the specific experimental situation 
at hand. In the standard text-book selective Von Neumann-Luders measurement (see e.g. 
|19p. the measured quantity is represented by a Hermitian operator S and the result of the 
measurement is given by an eigenvalue of S with probability depending on the current state. 
However several different scenarios and mathematical models of quantum measurements 
can be considered in different situations (see e.g. Therefore different definitions of 

observability may be appropriate and of physical interest in different cases. Nevertheless, 
there are several reasons to study observability for quantum mechanical control systems. 
From the viewpoint of the fundamental development of the theory, observability is one of the 
main concepts to be extended to quantum systems. It is related to the notion of input-output 
equivalence and therefore to the general question of modeling time varying Hamiltonians 2 . 
The problem of determining the state from the observation of a quorum of observables is 
an important one in quantum mechanics |16| . Techniques to find a set of observables which 
would determine the state without ambiguity have been extensively studied in quantum 
physics (see e.g. 0], |10j). Observability of quantum systems is also particularly important 
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2 Two models are input output equivalent if they produce the same output function for any input. Two 
input-output equivalent models cannot be distinguished by applying control inputs and observing the output 
and therefore modeling via input-output experiments may only be made up to equivalence classes of input- 
output equivalent models. This question is explored for networks of particles with spin in 0, 
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in view of the recent interest in implementing feedback at the quantum level (see e.g. [5], 
JBj, IZUi A feedback controller uses the knowledge on the current state to update 

the value of the control, i.e. it is of the form u = u(t,x). The knowledge of the state is 
obtained through the output and therefore an a priori knowledge of the extent to which 
information on the state can be obtained from the output is essential in the design of state 
feedback control scheme. 

In a recent paper [7|, a study was presented on the observability properties of quantum 
systems subject to nonselective measurement i.e. a measurement where either the result 
is not read or it is given by the expectation value of a given observable. The latter case 
is of interest in several experimental scenarios such as nuclear magnetic resonance where 
the output signal is averaged over a large number of quantum systems. In these cases, the 
definition and treatment of observability is simplified by the fact that one does not have to 
consider probabilities explicitly and natural definitions of observability can be given. In this 
paper we expand upon the treatment of |Z| for general measurements. A unified treatment 
for the various types of measurements is presented using notions of generalized measurement 
theory |Jj]. 

We shall be interested in the dynamics of finite dimensional quantum systems whose state 
is described by a density matrix p. We shall consider measurements occurring at discrete 
instants of time. In between two measurements, the evolution of p is governed by Liouville's 
equation (see e.g. |TU)) 

ip=[H(u(t)),p], (3) 

where the Hamiltonian H explicitly depends on a control u — u(t). In general, for nonselec- 
tive measurement the result can be assumed to be a linear function of the current state p. 
This is the case when one performs a Von Neumann-Luders measurement of the expectation 
value of a given observable S in which case the output y associated to a system © is given 

by 

V = Tr(Sp). (4) 

Another example is the indirect measurement discussed in detail in Section^ We shall treat 
the nonselective case in greater detail and then present some extensions to the selective case 
in Section H3 

The effect of nonselective measurements on the state p of the system can be described 
in general using the formalism of operations 0, [15]- I n particular, if M. is a measurable 
set of possible outcomes, upon measurement the state p is modified as 

p -» T(p) := \ <$> m {p)dm, (5) 

JM 

or 

p - T(p) := Y, *"»(p)> (6) 

according to whether M. is a continuous or discrete set respectively. The super-operators $ rn 
are called operations and, according to Kraus representation theorem |15| , can be expressed 
as 

k 

for a countable set of operators fl m k- 
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The plan of the paper is as follows. In Section |3] we give the basic definitions and results 
concerning the observability of quantum systems, under general nonselective measurement. 
The main concepts and results given in this section were already presented in |7j, however 
they are summarized here in a more compact form by introducing effective observables. In 
Section[3]we present some results for the special but important case of indirect measurement. 
These include an expression for the effective observable in this case and the derivation of 
the optimal measurement in terms of minimal disturbance on the state. This raises the 
question of whether there is a conflict between observability and minimal disturbance. In 
Section^] by a simple physical example we show this is not the case: in general we can have 
observability with a low disturbance of the system. In Section[3]the design of quantum state 
reconstruction is discussed and related to observability. Section presents an extension of 
the formalism to the case of selective measurement. 

2 Observability under general nonselective measurement 

If the output y of system Q is a linear function of the current state, as we assume here, it 
is always possible to express y as 



for some Hermitian matrix S e ff, which represents an effective observable. Without loss of 
generality, we can assume that S e ff has zero trace since a trace different from zero would 
only introduce a constant shift in the value of the output which does not play any role in 
our treatment. Alternatively, we could quotient all the subspaces (the observability spaces 
defined in (11 1)11 below) by span{il}. 

Denote by pk(t, u, p) the solution of © with initial condition p, control u at time t after 
k — 1 measurements, where, at every measurement, the state is modified as in l|H | -l(7 )l . Then, 
two states p\ and pi are called indistinguishable in k steps (or after k measurements) if, for 
every control u and time t 



A system is called observable in k steps if indistinguishability in k steps of p\ and p2 implies 
P\ = Pi- A system is called observable if it is observable in k steps for some k. 

As in the study of controllability (cf. |17|. [lip the dynamical Lie algebra associated 
to the quantum system @ plays a prominent role. The dynamical Lie algebra C is defined 
as the Lie algebra generated by span ue u{— iH(u)}, where U is the set of possible values for 
the control u. In order to express the conditions for observability in an arbitrary number 
of steps, under general nonselective measurement, we associate to the super-operator T a 
dual super-operator T* acting on observables S and defined from the requirement that, for 
every S and p, Tr{T* (S)p) — Tr(SJ-(p)). Then, we define generalized observability spaces 
Vfe, k = 0, 1, recursively as 



y(t)=Tr(S effP (t)), 



(8) 



Tr(S e ffPk(t,u,pi)) = Tr(S e ffpk(t,u,p 2 )). 



(9) 



V := span{iS e f f }, 



Vi:=0 



OO 



ad 3 c V , 



(10) 



V k := 



OO 



o4^*(v k _i), 



3=0 



where ad J c V is defined as spanned by all the repeated Lie brackets [Ri, [Hq, . . . , [Rj, iA 

and the Lie bracket is taken j times, R\, . . . , Rj € C and iA £ V. 

With these definitions, the main results of can be summarized as follows. 
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Theorem 1 System J3J) with output y in is observable in k steps if and only if 



V fe = su{n). (11) 

More in general, write p = pi + p% where p\ is the component of p in iVk 3 and p2 is the 
component along iVjjr where is the orthogonal complement ofVk in u(n). Then, we have 
the following decomposition of the dynamics 

pi = 

(12) 

p 2 = -i[H(u),p 2 ], 

and we have 

y(t) := Tr(S eff p(t)) = Tr(S effPl (t)). (13) 
Initial states are indistinguishable in k steps if and only if they differ by an element in iV^ . 

In several interesting scenarios, the measurement scheme has a 'repetition property' which 
can be defined by imposing that the operators Q m k in satisfy f2 m fcfi r ; = 6 mr dki^i m k, 
Vm,r £ M. and Vfc,L In these cases 

*m($ m (p)) = $ m(p) Vp, T 2 = T, and T* 2 = T*. 
Physically this means that a second measurement does not modify the state more than the 
first one. In these cases, it is easy to show that 

V fc -i C V k (14) 

so that states that are indistinguishable in k steps are also indistinguishable in k — 1 steps 
4 Moreover, because of the assumption of finite dimensionality, there exists a k such that 
Vfe = V^. for all k > k. An example is the standard Von Neumann-Luders measurement of 
the observable S. In this case S e ff = S. Expressing S as 

s = X>n,, (17) 

3 

where the Aj's are the eigenvalues of S and IL; are the orthogonal projections onto the 
corresponding eigenspaces which play the role of f2 m fc's. T is given by 

3 

In order to use the results of Theorem ^ we need to find an expression for T and S e f / 
which describe the particular measurement considered. In the following section we treat in 
detail the case of indirect measurement. 



3 vector space of Hermitian matrices obtained by multiplying by i the skew-Hermitian matrices in Vfe 
4 The proof uses an expression of S e ff in terms of effects F m defined in Section |H] When the output is 



an expectation value, then 



S e ff = ^ mFm - ( 15 ) 



Moreover using the expression for the effects 



^mfc^mfc, (16) 



and the repetition property, one has T* (S^f f) = S e f f and therefore Vi = T* (Vo). Vrj Q Vi and by induction 
one obtains 1141 . 
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3 Observability under indirect nonselective measurement 

In indirect measurement, the system evolves as in Q until it is in a state p$ and it is put 
in contact with a probe system whose initial state we denote by pp. The total system of 
system and probe at the beginning of the measurement process is in the state 

Ptot ■= Ps® Pp- (19) 

During the measurement process, of duration r, the total system evolves according to an 
Hamiltonian 

Htot ■= H{u)®l + g{t)A® B + 1® Hp. (20) 

The term Hp describes the dynamics of the probe system alone. The term g(t)A (g) B gives 
the interaction between probe and system, where g(t) is nonzero only during the interval 
[0,t]; 1 is the identity operator. It is usually assumed that, when the interaction is active, 
it represents the dominant term in the Hamiltonian Htot- Therefore we shall first assume 

Htot := g(t)A ® B. (21) 

At the end of the interval [0,r], an observable S is measured on the probe system, or 
equivalently an observable \ (&S is measured on the total system. In the following proposition 
we calculate an expression for S e f / with the Hamiltonian (|21|l . 

Proposition 3.1 With the above definitions and notations, for indirect measurement 

S eff = Y J A k Tr((ad k _ iB pp)s)—, (22) 

fe=0 

where 

G := / g{t)dt. (23) 



Proof. The solution of © with initial condition ptot in (|19f) and Hamiltonian Htot in 
(I21|l can be written at time r as 

PTOT ( T ) = e- iGA ® B Ps ® p P e* GA ® B . (24) 

Expanding, using the Campbell-Baker-Hausdorff formula, this can be written as 

°° G k 
Ptot(t) = ad A®-iB Ps®Pp-tj- (25) 

fc=0 

Now, it is easily seen by induction on k that every operator ad k i ^_ iB ps® Pp can be written 
in the form 

ad%_ iB PS ®PP = Y1 F i ® L J > ( 26 ) 

where the Fj 's are all operators of the form 

F r .= A k - l Ps A\ (27) 

for some I, < I < k, and 

2 fc 

Y l L j = ad k _ iB p P . (28) 
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Using i12.")|i and (|26|l in the expression of the output y, we obtain 

y = Tr(l <g> Sptot{t)) = 

00 G k 
Tr(l®sY,ad k m _ tBPs ®pp-^)= (29) 

fc=0 

00 „ fc 2 k 

EirE T ^^), 

fc=0 ' j=l 

where Fj and L 3 are defined in (|26|l . Using (|27|l . 128() and elementary properties of the trace 
we obtain 

00 2 k 

y = E -k\ Tr{Akps) E T MSLi) = (30) 

fe=0 ' j=l 

r / 00 f' k \ i 

Tr[(J2A k Tr P (ad k L lB p P S) — )p s . 

k=0 

As this has to hold for every p$, S e fj, in. y = Tr(S e ffps), is given by l|22|). □ 
We notice some features of the expression of S e ff l|22|l . 

Remark 3.2 Assume we retain only the terms up to first order in G. This is reasonable if 
the interaction is very quick and of small magnitude. Then we have 

S eff *Trp(ppS)l+Trp([-iB,pp]S)GA, (31) 

so that, if Tr{[—iB , p P ]S) ^ there is a one to one correspondence, in first approximation, 
between the values of the output and the value of the observable A, and therefore we can 
say that we are measuring A indirectly. 

Remark 3.3 In the special case where S and B are canonically conjugate observables on 
the probe, i.e. 

[B,S]=i^l, (32) 

with 76 R, the above correspondence between mean values of S e ff and A is exact. This 
is the case treated in @. In order to see this, consider the expression of S e ff (|22|l . From 
the property 

Tr P {{ad k L lB p P )S) = (-l) k Tr((ad k _ tB S) PP ) (33) 

and Q32|l ■ we obtain that all the terms in the sum corresponding to fc > 2 are zero as ad k _ iB S 
is zero in these cases. Therefore <5*e// reduces to 

S eff = Trp{p P S)\ + jGA, (34) 

where 7 is the one in 1)32(1 and we have used the fact that Tr(p P ) = 1. 

Remark 3.4 In some cases, it is not appropriate to neglect the term containing H(u) in 
(|20|l . In these cases, it is not possible, in general, to obtain a simple expression of S e ff as 
in (|22|l . However Remark 13.21 above still holds true, assuming the g{t) is a simple square 
function in [0, r] so that G — r, and u is constant in [0,r]. Notice that we can write, 
generalizing J2U), 

PTOt(t) = ad -tH(u)®l+A®-tB PS ® PP-jj - (35) 
k=0 



G 



and expanding Ptot{t) to first order in r, we have 

Ptot ~ PS® Pp + ad-iH(u)m+A®-iB Ps <8> Pp t. (36) 

Calculating Tr(l <g> Spror), using the fact that [— iH(u) ® 1, ps (8 pp] = [— iH(u), ps] ® pp, 
and that Tr([— iH(u), ps}) — 0, we obtain the same expression for S'e// as in i|31|) . with G 
replaced by r. 

Remark 3.5 The expression of S e ff does not depend on the probe being finite dimensional. 

The action on ps (J 7 in (J5J, JJ]J) after an indirect measurement is given by 

Ps - Hps) := Tr P (e A ®~ WG p s ® pp e A ^ SG ), (37) 

and it is independent of the observable S measured. This is easily seen by using the properties 
of the partial trace Trp and the expression of S <|17|) in terms of projections. We refer to |H] 
(Section 2.4.6) for a discussion of how this operation on the state can be rewritten according 
to Kraus representation theorem as in ©-10. 

In J22J), there is a dependence of S e ff on the initial state of the probe. As a consequence, 
it could be possible to modify the observability property for the system by suitably choosing 
pp. However, the disturbance induced on the system depends on pp as well, and it is inter- 
esting to investigate whether there is a conflict between observability and low disturbance 
of the system. 

With this motivation in mind, we provide here an analysis of the disturbance on the state 
while performing an indirect non selective measurement and show how to find the initial 
state of the probe which gives the (worst case) minimal disturbance. Using this result, we 
shall show in the next section, with an example, that there is in general no conflict between 
observability and minimal disturbance. 

We consider, as a measure of the disturbance on the state ps, the trace norm 

d := \\F(Ps) - Ps(0)\\ = [Tr(T(p s )-p s ) 2 }K (38) 

expressing the distance between the initial state ps and the final one, J-(ps)- If we fix all 
the parameters of the measurement process, the disturbance d will in general be a convex 
function of ps- Since ps varies on a convex and compact set, the set of all the density 
matrices, the maximum will in general be achieved on the boundary i.e. it will be a pure 
state. We shall now show how it is possible to find this worst case pure state in the small 
time approximation in the case where all the terms in 1121 III are possibly different from zero 
(and u is constant). After that, we will derive the corresponding distance d, depending on 
pp. Then, it will be immediate to find the initial state of the probe which gives the minimum 
for d. In the above situation, neglecting higher order terms in r, d 2 can be written as 

d 2 = -t 2 Tr([H(u) + Trp(Bp P )A, p s (0)]) 2 ■ (39) 

If we set 

X :=H(u)+Trp(B PP )A, (40) 

we can write d 2 as 

d 2 = 2T 2 Tr(X 2 p 2 - XpsXps), (41) 

where we write ps for ps(0) as there is no possibility of confusion. As an orthonormal basis 
for the Hilbert space of the system, we choose the eigenvectors of the Hermitian operator 
X, \(j>k), k = 1, ...,rc., so 

n 

^ = 5> fc |<M(0 fc | (42) 
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and Xk are the real eigenvalues of X. Since the worst case ps is a pure state, we can write 
ps = f° r some \ip) = Efc^l^fc) wnere the n coefficients Tk completely specify ps- 

They can be assumed real by suitably redefining the eigenvectors \<fik)- We have the further 
constraint r\ — 1 since Trps = 1. To determine the worst case ps, we rewrite (|41 1) as a 
function of the coefficients 

n n 

S = 2r^ X lrl-Y^^i{r k r l f) (43) 

k=l kl=l 

and, rearranging the terms, 

d 2 =2r 2 (^{x k -x 3 ) 2 {r k r ] ) 2 ). (44) 

k>j 

We now maximize d 2 with respect to the n parameters r k using the Lagrange method: 

d r J(ri,...,r n ) = forZ = l,...,n 

(45) 

dxf(n, . . . ,r„) = 0, 
where 

n 

/(n,..., r„)=d 2 + A(^r fc -l) (46) 
fe=l 

and A is the Lagrange multiplier. More explicitly, 

n [12jjti( x i ~ x jT r ) + A) = for / = 1, . . . ,n 

(47) 

System (|47|l always admits a solution since the function d 2 is continuous over the compact 
set of pure density matrices. In the next section we will explicitly compute ps in a particular 
case. Without solving (|47|) in the general case, we summarize our discussion in the following 
theorem. 

Theorem 2 The worst case disturbance in a small time approximation is given by d 2 in 
\44\l ; where (ri, . . . , r n ) are the solution of system Therefore given u, A and B in the 

definition of X , the initial state of the probe which minimizes the worst case error has to be 
chosen so as to minimize this d 2 . 

4 Observability and minimal disturbance 

As a concrete example of observability under an indirect measurement, we consider the 
simple case of two-dimensional system and probe. The system is a qubit with external 
control u affecting a two-components magnetic field, for example 

H{u) = E x (u)a x + E y (u)<j y . (48) 

We assume a piecewise constant control u £ {ui, U2] that flips the magnetic field directions 
x and y, that is E x (u{) — E, E y {u\) — and E x {u2) — 0, E y {u-i) — E. We use a second 
qubit as probe and we let it interact with the system for a short time r in which the free 
evolution (|48|l can be neglected. To get information about the initial state ps we measure 

5 = a z on the probe. Assuming a simple Ising model of interaction, A — a y and B = a x , 
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the effective observable S e tf can be explicitly computed. Splitting the sum in llL'l'l) in even 
and odd indices, using (|3*3^l and considering that 

L for fc even, 
A K = i (49) 



a y for k odd 



and 



we find that 



{(— l) fe / 2 2 k a z for k even, 

(50) 
(-l) 1+k / 2 2 k+1 cr y for fc odd 

S eff =Trp(cr zP p) cos 2G 1 + Tr P (a y p P ) sax{2G)a y . (51) 



Remark 4.1 The observability properties of our system strongly depend on the initial state 
of the probe pp. Suppose that Tr{a y pp) = 0; in such a case S e ff = and the observability 
spaces Vk contain only the null vector. Then, for any k the system is not observable and 
the states are all indistinguishable. On the other hand, suppose Tr(a y pp) ^ 0. In such a 
case S e ff — Trioypp) sm2Go y and Vk — su(2) for all fc, and the system is observable in fc 
steps, for every fc. 

In some cases, it is not appropriate to neglect the free evolution of the system. However, 
following Remark 13.21 we can explicitly evaluate the effective observable: 

S e ff = Tr P (a z pp) 1 + 2TTr P (a y p P ) a y (52) 

where r is the time of interaction (assumed small) and g(t) is a square function. Remark 
14. II holds true in that case as well. 

We now determine the minimal disturbing probe described in Theorem [21 We assumed 
that during the time interval r the control does not change, and its actual value is relevant 
in order to find the minimal disturbing probe. In our example, (14711 becomes 

ri({x 2 - Xl ) 2 rl + \) =0 

r 2 ((x 2 -x 1 ) 2 r 2 + A) =0 (53) 

r l + r 2 = 1- 

where x\, x 2 are the eigenvectors of X and they depend on u. Solving (|53|l we find the worst 
case ps- 

ps = l (l^i)^il + \h)(<f>2\ ± \4>i)(<h\ ± \(h)(M (54) 

leading to d 2 = (x 2 — £i) 2 /4. For u = u u x 2 — X\ — 2[E 2 + (Tr P (a x p P )) 2 ], for u = u 2 , 
x 2 — X\ = 2(E + Trp(a x pp)). Then, the minimally disturbing probe must satisfy 

Tr P (a x p P ) = foru = ui, 

(55) 

Trp(tj x pp) = max{-i?, —1} for u = u 2 . 

In both cases there is not a conflict between observability and minimal disturbance (see 
Remark |4~T|) . 
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5 Observability and state reconstruction 



We present in this section a system theoretic treatment of the problem of state determination 
for the system (|3J) with output Q. In systems and control theory, for a continuous time 
system such as @, under observability conditions, the (initial) state is determined from a 
continuous reading of the output. From a physics point of view, a continuous monitoring of 
the output will introduce a back action on the state of the quantum system and therefore it 
will render invalid the model ©. However, this scheme is of interest for quantum systems 
in situations like the following. Assume we want to determine the unknown (initial) state 
and we have many copies of the same system. We perform a nonselective measurement on 
each copy at slightly different times so as to simulate a continuous measurement. The data 
so obtained can then be used by the observer to reconstruct the state of the system (without 
measurement back-action). 

With this motivation in mind, a method for reconstructing the initial state can be ob- 
tained by adapting to our case techniques for time varying linear systems |14j . Observability 
(in one step) is a necessary and sufficient condition for reconstructing the initial state from 
a reading of the output. In fact, if the system is not observable, then it is not possible 
to discern between two indistinguishable initial states. Viceversa, assume the system is 
observable. Then, we have that [7] 

{X*iSX\X G e £ } = su(n). (56) 

This means that we can choose a control u, so that, for the corresponding solution X u of 
Schrodinger operator equation 

X = -iH(u)X, X(0) = I, (57) 

the n 2 — 1 elements of the matrix X*SX U (namely the real functions composing the matrix 
modulo the fact that this matrix is Hermitian) are linearly independent. e c is the Lie group 
of all the matrices X m for which there exists a control steering X in (|57(l from the identity 
to X m . We can select n 2 — 1 matrices Xi 7 ...,X n 2_ 1 so that X* 1 S l Xi,...,X* 2 _ 1 S'X n 2_ 1 are 
linearly independent and then concatenate the controls steering the matrix X in (|57|l to 
Xi, X 2 X*, X 3 X|,...,X n 2_ 1 X* 2 _ 2 . Now assume that, in the control interval [0,T], the 
(significant) real entries of X*SX U are linearly independent and define the linear operator 
W which maps nx n Hermitian matrices with zero trace into nx n Hermitian matrices with 
zero trace as follows 

W u (p ) := / X:(t)SX u (t)Tr(x:(t)SX u (t)p )dt. (58) 
Jo 

The operator W u has the following property. 

Proposition 5.1 If the n 2 — 1 real functions composing X*SX U are linearly independent 
then W u has rank n 2 — 1 and therefore it has an inverse W" 1 

Proof. This follow from the well known fact that (see e.g. Section 9.2.1) m functions 
lj = lj(t), j = 1, m are linearly independent in an interval [0, T] if and only if the matrix 

9ij := f li{t)lj(t)dt (59) 
Jo 

is nonsingular. In our case, if we order the n 2 — 1 elements of po by row and then with real 
and imaginary part i.e. as p 0l l , Re(p' 0l 2 ), Im(p 0l 2 ),...,i?e(p 0( „_ 1) n ), /m(p ( n -i),n)) Po n ,n 
and the entries of X*SX U in the same way, the matrix which represents the linear application 
W u has the form (|5*97l where k are the elements of X*SX U and therefore it is invertible. □ 
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Now, from formula we obtain 



y(t) = Tr(X*SX u (p - -J nxn )), (60) 

71 



and therefore 



X*{t)SX u (t)y{t)dt = / X^(t)5X u (t)Tr(X*5X„(p --i"„xn))rft. (61) 



Therefore, using the definition of W u (|58|l . we have the following formula for the reconstruc- 
tion of the initial state po, 

Po = i/„ x „ + W- 1 (jf X*Jt)SX u (t)d?j . (62) 

Formula (|62|l represents a system theoretic alternative to methods for quantum state to- 
mography. We summarize the discussion in the following theorem. 

Theorem 3 Consider system 0) with output Q). If the system is observable (in one step), 
then there exists a control such that formula ffi^J) gives the initial state. 

An alternative to the 'static' state reconstruction formula Ij62(l is the design of an asymp- 
totic observer namely a dynamical system which uses only a reading of the output and whose 
state asymptotically converges to the actual state of the system. We present in the rest of 
this section proposal for such an asymptotic observer which is inspired the treatment for 
linear time varying systems in |12| . We notice that for static state reconstruction we im- 
posed a requirement on the control which impled that the operator defined in formula Ij58(l 
has full rank (cf. Proposition l5.1fl . For an asymptotic observer which estimates the state as 
t — > oo, we need to impose that this property is somehow uniform for every t as t — > oo. To 
make this more precise we define a time dependent, linear symmetric operator on Hermitian 
matrices Vt as follows (we omit for notational simplicity the dependence on the control u). 
Let U := U{t) be the solution of Schrodinger operator equation 

U = iH(u)U, 17(0) = I, (63) 

where H{u) is the same as in (0) and (|57[1 . We define 

V t (A) :=V^' a (A) := f e- M{t - T) U{t)U*{T)SU{T)U*{t)Tr(U{t)U*{T)SU{T)U*{t)A)d T , 

(64) 

with a > and M > 0. We assume that the control u is such that there exists a a > 0, 
such that, for every t > a 



ai Tr{A 2 )< f (Tr(U(t)U*(T)SU(T)U*(t)A)) 2 dr 

J t-a 



< a 2 Tr(A 2 ), (65) 



for some positive constants a\ and ai. 

We choose the same a in the definition (|64fl and our assumption implies that 

ai e- Ma Tr(A 2 ) < Tr(AV t (A)) < a 2 Tr(A 2 ), (66) 

and that Vt is nonsingular, so that we can define the inverse operator P t _ . Moreover from 
the definition Q64f) . we obtain (A constant) 



^P t (A) = STr(SA)-e- Ma U(t)U*(t^a)SU(t-<7)U*{t)Tr(U(t)U*(t-a)SU(t-a)U*(t)A) 

(67) 
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-MV t (A) + [iH(u),V t (A)]+V t {[iH{u),A]). 
We consider the state observer for system © 

j t P=[-iH{u),p] - hpt\S)Tr{S{p-p)) = [-iH(u),p] \r^{S)(y Tr(Sp)). (68) 

Here p is the estimate of the actual state p and we are going to show that A(t) :— p(t) — p(f) 
tends to zero as t — > oo. By subtracting l|68[l from Q, we obtain the differential equation 
for A, 

A = [-iH(u), A] - lv t (S)Tr(SA). (69) 

Define the Lyapunov candidate function V := V(t, A) = Tr{APt{A)) 1 which according 
to l|6*6"|) satisfies 

a[Tr{A 2 ) < V(t, A) < a' 2 Tr(A 2 ), (70) 

for appropriate positive constants a[ and a' 2 . Moreover we can calculate -^V(t, A(t)). By 
using (Hp and $7$ along with the property Tr(P t (A)'P i r 1 (S)) = Tr(AS), we obtain 

j t V(t,A)<-MV(t,A), (71) 

and therefore it follows from Lyapunov second method |!5] that system (|69|l is asymptotically 
stable, and therefore A tends to zero as t — > oo. Notice that our Lyapunov is only defined 
for t sufficiently large (t > a) however this does not change the stability analysis as all 
the quantities considered are guaranteed to be bounded over a finite interval of time. We 
conclude with the following Theorem. 

Theorem 4 Consider system with output Assume that the control u satisfies the 
condition \65}) . Then system f6'^|) is an asymptotic observer for |3JJ. 



6 Some extensions to selective measurement 

In this section, we discuss how the theory described above for nonselective measurement 
extends to selective measurement. There is no difficulty in doing this in the most general 
case namely in the context of the generalized measurement theory of operations and effects 
[Hj. According to this theory, given a measurement scheme, to every result m is associated a 
positive operator F m , called an effect. If p is the current state of the system, the probability 
of obtaining the result m (or of an event m to occur) is 

P(m) = Tr(F m p). (72) 

After a result m (or, more generally an event m) has occurred, the state is modified according 
to 

p^Pim)- 1 ^), (73) 

where the positive super-operators $ m are the same operations as in Q) and Tr{^ m (p)) — 
P(m) = Tr(F m p). Two initial states p\ and p 2 are said to be indistinguishable in k steps, 
in selective measurement, if they give every possible result with the same probability at the 
fc— th measurement, for every choice of the control u. In formulas (cf. ©) 

Tr(F mPk (t,u,p 1 ))=Tr(F m p k (t,u,p 2 )) Vm e M, (74) 

where M. is the set of possible results (events). Let Pk(m) be the probability of having 
the result m at the k— th measurement and let P(mi, ...,mfc) be the joint probability of 
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having result mi at the first step, mi at the second step and so on. Also, indicate by 
P k (m k \mi, . . .mk-i) the conditional probability of having at the k— th measurement, 
given mi,. . . m k as ordered results of the previous measurements. By use of the formula 

P k (m)= ^2 Pk{rn\mi,...,mk-i)P{rni,...,m k -i), (75) 

i»l...nin 

and repeated use of Bayes' formula 

P(mi, ...,to*_i) = (76) 

P fc _i(m fc _i|mi, ...,m k _ 2 )P(m 1 , ...,m fe _ 2 ), 
we can write (m) starting from an initial condition p as 

P fe (m)= ]T rr(F„,Xjfe(4 niJk _ 1 (Jr fc _i 

mi...mj._i 

. . . (* mi (*iA>*i)) • . (77) 

where Xj, j = 1, . . . , k is the evolution solution of the Schrodinger operator equation l|57|) 
in the interval between the (j — l)-th measurement and the j— th measurement. Using <|77|) 
and using the linearity of the operators <!>,„, we can rewrite Pk(m) as 

P fe (m) = Tr(F m X k T{X k _ 1 ...F{X lPo Xt)...X* k _ l )Xt) (78) 

where T is defined in © . From this point on the theory goes as in [7] and the result is an 
extension of Theorem ^ In particular, one defines the 'selective' observability spaces (cf. 

:= span meM {iF m }, Vf ! := ®~ ad 3 c {Vf), 

(79) 

vr ! :=e^o«4^*(v^i), 

and Theorem ^ extends by replacing nonselective observability with selective observability 
and the spaces V with the spaces V sel 5 . 

The remarks following Theorem ^ on the implications of the repetition property also 
extend with only minor formal modifications. In the particular case of the standard Von 
Neumann-Luders measurement, the observable S is written in terms of the projectors 
and the eigenvalues A as 

s - xn ^ ( 8 °) 

and the above theory holds with H\ playing the role of the effects F m . 

Remark 6.1 The observability space Vq s1 does, in general, include the observability space 
Vo and therefore the same is true for the observability spaces V k el and Vk ■ This implies that 
nonselective observability implies selective observability, as it is intuitive but not viceversa. 
Consider as a specific example a spin 1/2 particle for which the z component of the spin is 
measured, with a Von Neumann-Luders measurement. In this case Vo = span {ia z }, while 
Vq c1 = span {i<Jz, il}, with a z the z-Pauli matrix. So up to the span of il the two subspaccs 
are the same and the observability properties in the selective and non selective case are the 
same and depend on the dynamics. However consider a spin 1 on which we perform a Von 

5 Condition 1161 of Theorem^needs to be slightly modified as the effects F m do not necessarily have zero 
trace, by replacing Vfe with Vk/span{il} or by making all the effects traceless. 
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Neumann-Liiders measurement of the spin along the z direction. In this case, using the 
representations of the spin angular momentum calculated for a spin 1 (see e.g. ^5] Section 
3.5), we have 

(i \ 

Vo = span{ \ }, (81) 
\0 -i) 

and 

( i 0\ /0 \ 
Vo e Vspan{«l} = span{ -i , i }. (82) 
\0 0/ \0 -ij 

Therefore Vo and Vq £1 / span{il} do not coincide in this case. In particular, if we consider 
the dynamics determined by a time- varying control electro-magnetic field in the x — y plane, 
we have that Vi is spanned by the three dimensional Lie algebra representation of su(2), 
while Vq c1 I ' span{il} is equal to su(3), as one can easily verify by repeated Lie brackets. 
Therefore, in this case, we have selective observability but not non selective observability, 
in one step. The situation is the same if we consider observability in k steps as Vi = Vk for 
every k. 



7 Conclusions 

This paper has presented a collection of results on the observability of quantum systems 
with emphasis on the case of nonselective measurement. In particular 

1. Using the formalism of generalized measurement and of effects and operations we 
have extended the basic definitions and criteria of observability to the case of general 
measurement by introducing an effective observable. 

2. We have derived a general expression for the effective observable in the case of Von 
Neumann indirect measurement. 

3. In the case of indirect measurement, we have derived an expression for the state of the 
probe which would introduce the minimum disturbance in the state to be measured. 
We have showed that the requirement of a minimal disturbing probe does not in 
general compromise the observability properties of the resulting system and therefore 
the amount of information obtained on the state by the measurement of the output. 

4. We have presented two system theoretic methods to reconstruct the state by a measure- 
ment of the expectation value of an appropriate observable. One of them is through an 
integral formula and uses readings over a finite interval of time. The other is through 
an asymptotic observer whose state converges to the state of the measured system. 

5. We have extended the basic definitions and observability criteria to selective measure- 
ments. 

We believe that the system theoretic approach to quantum state determination is worth 
being further investigated. Extensions of our definitions and results to continuous measure- 
ments, optimization of the methods for state determination in specific settings, applications 
of observer design in closed loop quantum systems are only few possible subjects for future 
research. 
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